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We study the Za\ decay mode of a heavy CP-even Higgs boson of the NMSSM, hi, where a\ is 
the lightest CP-odd Higgs state of this scenario, the former produced in association with a bottom- 
antibottom pair, and find that, despite small event rates, a significant (in fact essentially background 
free) signal should be extractable at the LHC with very high luminosity, so long that a Z — \ jj 
(where j represents a jet) and ai — > t + t~ final state is exploited, in presence of 6-tagging. 

PACS numbers: 



The Next-to- Minimal Supersymmetric Standard 
Model (NMSSM) Q, owing to the introduction of an 
extra complex Higgs singlet field, which only couples 
to the two MSSM-type Higgs doublets, embeds a 
Higgs sector which comprises a total of seven mass 
eigenstates: a charged pair h , three CP-even Higgses 
^1,2,3 {mhi < m h 2 < m h a ) and two CP-odd Higgses ai^ 
(m 0l < m 02 ). 

Consequently, Higgs phenomenology in the NMSSM 
can be very different from that of the MSSM. As a key 
example, over the past few years, there have been several 
attempts to extend the so-called 'No-lose theorem' of the 
MSSM - stating that at least one MSSM Higgs boson 
should be observed through the usual SM-like produc- 
tion and decay channels at the Large Hadron Collider 
(LHC) throughout the entire MSSM parameter space @ 
- to the case of the NMSSM From this perspec- 

tive, it was realised that at least one NMSSM Higgs boson 
should remain observable at the LHC over the NMSSM 
parameter space that does not allow any Higgs-to-Higgs 
decay mode. In contrast, when a light non-singlet (and, 
therefore, potentially visible) CP-even Higgs boson, hi 
or h 2 , decays mainly to two very light CP-odd Higgs 
bosons, - ► o-io-i, one may not have any Higgs sig- 
nal of statistical significance at the LHC [6j. In fact, 
further violations to the theorem may well occur if one 
enables Higgs-to-sparticle decays (e.g., into neutralino 
pairs, yielding invisible Higgs signals) [23j |. While there 
is no conclusive evidence on whether a 'No-lose theorem' 
can be proved for the NMSSM at the LHC, there has also 
been put forward an orthogonal approach: to see when 
a, so to say, 'More-to-gain theorem' for the LHC 0-Q 
can be formulated within the NMSSM. That is, whether 
there exist regions of the NMSSM parameter space where 
more and/or different Higgs states are visible at the LHC 
than those available within the MSSM. 

In our attempt to overview both such possibilities, we 
assume here a light CP-odd Higgs boson, m, emerging 
from the decay h 2 — > Za\, where the heavy CP-even 
state, h 2 , is produced in association with 6-quark pairs 
at the LHC. This work complements the one carried out 
in a previous paper where the light a\ state is produced 
in pairs from decays of not only h 2 states but also hi 
ones, always from the aforementioned production mode 



[1 01 ] . Direct a\ production in association with bb pairs 
was considered in [ll[ for the case of a\ — > t+t~ as well 
as a\ — > 77 decays, in [l2| for the case of the a\ — > 
channel and in |13| for the a\ — > bb mode. 

For our study of the NMSSM Higgs sector we have 
used NMSSMTools 0, which is a numerical pack- 
age computing the masses, couplings and decay widths of 
all the Higgs bosons of the NMSSM in terms of its input 
parameters at the Electro- Weak (EW) scale. NMSSM- 
Tools also takes into account theoretical constraints as 
well as experimental limits, including the unconventional 
channels relevant for the NMSSM. 

Here, instead of postulating unification or taking into 
account the SUSY breaking mechanism, we fix the soft 
SUSY breaking terms to a very high value, so that they 
give a small or no contribution at all to the outputs of 
the parameter scans. Consequently, we are left with six 
free parameters at the EW scale, uniquely defining the 
NMSSM Higgs sector at tree-level. Our parameter space 
is in particular identified through the Yukawa couplings 
A and k, the soft trilinear terms A\ and A K , plus tan/3 
(the ratio of the Vacuum Expectation Values (VEVs) of 
the two Higgs doublets) and /i e s = X(S) (where (S) is 
the VEV of the Higgs singlet). 

In order to make a comprehensive study of the NMSSM 
parameter space, we have used NMSSMTools to scan over 
the aforementioned six parameters taken in the following 
intervals: 

0.0001 < A < 0.7, 0<k<0.65, 1.6 < tan/3 < 54, 
100 GeV < fi cS < 1 TeV, - 1 TeV < A\ < 1 TeV, 
-10 GeV <A K <0. 

Soft terms which are fixed in the scan include: 

• m Qa = m Ua = m Ds = m Ls = m Ea = I TeV, 

• A Us = A Ds = A E3 = 1.2 TeV, 

• m,Q — mr/ = tod = to^ = the = 1 TeV, 

• Mi = M 2 = M 3 = 1.5 TeV. 

We have finally performed our scan over 2xl0 7 ran- 
domly selected points in the specified parameter space. 
The points which violate the constraints (either theo- 
retical or experimental) are automatically eliminated by 
NMSSMTools. 
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The surviving data points are then used to determine 
the cross-sections for NMSSM Higgs hadro-production 
by using CalcHEP [H] for signals and MadGraph 
for backgrounds. As the SUSY mass scales have been 
arbitrarily set well above the EW one (see above), the 
NMSSM Higgs production modes exploitable in simula- 
tions at the LHC are those involving couplings to heavy 
ordinary matter only. Amongst the productions channels 
onset by the latter, we focus here on 



pp(qq,gg) -^bb h 2 , 



(1) 



i.e., Higgs production in association with a 6-quark pair, 
followed by 



(2) 



(incidentally, notice that h± is never heavy enough to 
enable sizable Za\ decays) and then 



Zai {jj){r + T ) (where j = jet). 



(3) 



Notice that we discard here the possibility of a\ — > 
decays, for two reasons: on the one hand, the mass re- 
gion below the t + t~ threshold is severely constrained 
(see and references therein); on the other hand, the 
decay rates are « (mP ole /mP° le ) 2 times suppressed 
with respect to the t + t~ ones, so that they would be 
numerically irrelevant (see forthcoming figures). As for 
a± — > bb decays, while being in turn more numerous 
than the t + t~ ones (above the 2mb threshold) by an 
amount oc (mb(mh 2 ) /mf olc ) 2 , they are overwhelmed by 
QCD backgrounds in Z +jet final states. Notice that the 
cti masses of relevance to this study will typically range 
between 15 and 60 GeV or so. 

We adopt herein CTEQ6L [H| as parton distribution 
functions, with scale Q = y/1, the centre-of-mass energy 
at parton level, for all processes computed. We finally 
assume y/s = 14 TeV throughout for the LHC energy. 
Also, in our numerical analyses, we have taken mb(mb) = 



4.214 GeV and 



pole 



171.4 GeV for the (running) 



bottom- and (pole) top-quark mass, respectively, while 
we have input m? T olc = 1.777 GeV and mj* = 0.1057 
GeV for the (pole) tau- and (pole) muon-lepton mass, 
respectively. 

As an initial step towards the analysis of the data, 
we have computed the production cross-section a(pp — > 
66/12) times the decay BK(h 2 — > Za\) against the BR 
itself, see Fig. 1. It is clear from this plot that most of 
the NMSSM parameter space yields rather small rates: 
in fact, only when the BR(h 2 — > Zai) is large the overall 
process (jTJ)— ([2]) can offer some chances of detection. In 
this scatter plot, one can notice a population of points 
with cross-section of up to several pb's. This region of 
parameter space is the only one exploitable at the LHC 
and reflects a well defined setup for the six input parame- 
ters of the NMSSM. Generally, such points correspond to 
generic and comfortably perturbative A and k (between 
0.05 and 0.15), very large tan/3 (above 35), rather large 




FIG. 1: a(pp ->■ bbh 2 ) BR(h 2 ->■ Zai) vs. BR{h 2 -¥ Zai) 
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FIG. 2: a{pp -> bbh 2 ) BR(h 2 -> Zai) BR(ai ->■ t+t~ 
BR(Z — >■ jj) vs rrih 2 . 



(and positive) /i e s and large (and negative) A\ (both 
between 200 and 900 GeV) and slightly negative A K (be- 
tween — 10 and —1 GeV). 

The rates in the previous figure refer to on-shell Z and 
ai though. One should clearly extract these from suitable 
decays, inevitably coming with a decay probability less 
than 1. Here, as intimated earlier on, we attempt the 
case of hadronic decays of the gauge boson and r-decays 
of the CP-odd Higgs boson. The corresponding signal 
yield for the ensuing final state is displayed in Fig. 2, 
now mapped against m^ 2 . After folding in the above 
decay probabilities, one is left with event rates at the 
0(100) fb level at the most. While clearly this number 
is not very large, signal events may still be detectable at 
planned LHC luminosities, especially if the background 
can be successfully reduced to manageable levels. Before 
proceeding to do so, it is worth mentioning that this kind 
of event rates are only found for h 2 masses well above 250 
GeV. 

We perform next a partonic signal-to-background 
(S/B) analysis. Amongst the backgrounds, we con- 



3 



Point 


A 


K 


tan/3 




1 


0.11784333E+00 


0.99759129E-01 


0.45413385E+02 


2 


0.36954734E-01 


0.74106016E-01 


0.35731787E+02 


3 


0.55822718E-01 


0.41921611E-01 


0.47269331E+02 


4 


0.15630654E+00 


0.84945744E-01 


0.48122165E+02 


Point 


Mcff [GcV] 


A x [GeV] 


A K [GeV] 


1 


0.65363068E+03 


-0.56108302E+03 


-0.97511471E+01 


2 


0.33852788E+03 


-0.68364824E+03 


-0.13135077E+01 


3 


0.38990121E+03 


-0.29909449E+03 


-0.53814631E+01 


4 


0.42293377E+03 


-0.23627361E+03 


-0.74575765E+01 


Point 


m at [GeV] 


m h2 [GeV] 


r oi [GeV] 


T h2 [GeV] 




1 


44.1 


272.4 


0.205 


6.57 


2 


15.9 


261.4 


0.01879 


4.972 


3 


61.5 


275.6 


0.00293 


7.86 


4 


31.9 


287.9 


0.0674 


8.37 






Point 


cr(pp bbhi) BR(h 2 -> Zai) [fb] 


BR(ai ->■ r+r") 




1 


1286.35 


0.141549576 


2 


377.3 


0.07615127 


3 


109.0 


0.125658544 


4 


447.1 


0.111203597 



bbh 2 -» bbZa : -> bbjjT+r- 



TABLE I: The NMSSM benchmark points used in the S/B 
analysis. 



sider here what we verified to be the dominant one, 
the irreducible noise induced by pp — > bbZr + r~ chan- 
nels (amounting to 60 Feynman diagrams in the unitary 
gauge). As benchmarks for the S/B analysis we have 
chosen four points in the above parameter space region 
(see Fig. 2), as illustrative examples. These benchmarks 
are given in Tab. I. 

In the light of the cross-section and decay rates in the 
table, and recalling that BR(Z — > jj) ps 70%, the total 
inclusive cross-section for the signal in (JTJ-([3]) for point 
1(2)[3]{4} is 127(20)[10]{35} fb. The SM irreducible 
background inclusive rate is 7.6 fb. After implementing 
the following acceptance cuts [24| 

AR(i,j) >0.4 = b,b,j,j,T + T~), 
|r?(i)|<2.5 (i = b,b,j,j,T + T-), 



Pr{i) > 15 GeV (i — b, b, j, j, r + , r ), 
and the selection ones as well 



(4) 



\M H -M Z \ < 15 GeV, 



\M T+T - 



i ai \ < 15 GeV (5) 



(wherein m ai need not be known beforehand, as it would 
be seen in the reconstructed M T + T ~ distribution: here 
it assumes the values of Tab. I), we obtain for the sig- 
nal cross-sections the values 2.27(0.12)[0.21]{0.65} fb for 
point 1(2) [3] {4} so that the typical signal efficiency isap- 
proximately 1.78(0.57) [2.21]{1.87}%, respectively (25|). 



Point 1: solid/black 

Point 2: dashed/red 

Point 3: dotted/blue 

Point 4: dot— dashed/green 



100 
pT(r) (GeV) 



FIG. 3: The (reconstructed) r transverse momentum distri- 
bution for the four benchmark points in Tab. I. 



The selection efficiency strongly depends on the m ai 
value, see Fig. 3, in the sense that the lighter a\ the 
softer its decay products (the r + r~ pair), so that the 
Pt{t + ,t~) selects fewer signal events: recall in fact that 
one has m ai = 44.1(15. 9)[61.5]{31. 9} GeV for point 
1(2) [3] {4}. Notice that after all such constraints the 
SM irreducible background is essentially removed alto- 
gether, as its rate goes down dramatically, to O(10~ 4 ) 
fb for all points 1-4. The distributions in invariant mass 
rrijj T + T - , that one could obtain after reconstructing the r 
decay products, for the four benchmark points, are given 
in Fig. 4, showing the Breit-Wigner peaks signaling the 
h-2 —> Zai resonances. Integrating over the entire spec- 
trum, one obtains 2274(115)[213]{653} signal events for 
point 1(2) [3] {4}, essentially background free, assuming 
1000 fb" 1 of integrated luminosity. 

In summary, we have proven that there exists a small 
but well defined region of parameter space where the hi 
and a\ states of the NMSSM, both with a mixed sin- 
glet and doublet nature, could potentially be detected at 
the LHC if 250 GeV < m h2 < 300 GeV and 15 GeV 
m ai ^ 60 GeV, in the h 2 — > Zai —> Ht + t~ mode, 
when the CP-even Higgs state is produced in associa- 
tion with a bb pair for rather large tan/3, large (and posi- 
tive) jU e ff; large (and negative) A\, and slightly negative 
A K , for typically perturbative values of A and k. After 
a realistic S/B analysis at parton level, we have in fact 
produced results showing that the extraction of such a 
signal above the dominant irreducible SM background 
should be feasible using standard reconstruction tech- 
niques [2(| HH and large LHC luminosities, i.e., after sev- 
eral years of running at design values or rather promptly 
at the Super-LHC [23]. While more refined analyses, in- 
corporating r-decays, parton shower, hadronisation and 
detector effects, are needed in order to delineate the true 
discovery potential of the LHC over the actual NMSSM 
parameter space, we are confident that our results are a 
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bbh s -> bbZa 1 -> bbjjr+T- 



- Point 1: solid/black 
Point 2: dashed/red 
Point 3: dotted/blue 

: Point 4: dot-dashed/green 
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FIG. 4: The (reconstructed) /12 mass peaks in terms of the 
jjr + r~ invariant mass distribution for the four benchmark 
points in Tab. I. We plot the number of events after an inte- 
grated luminosity L = 1000 fb _1 . 



step in the right direction to both: (i) prove the existence 
of a 'More-to-gain theorem' at the CERN collider for the 
NMSSM with respect to the MSSM (as Higgs -> Z Higgs' 
signals are only possible in the latter scenario in parame- 
ter space regions already excluded by experimental data) 
and (ii) to establish a 'No-lose theorem' for the NMSSM 
at the LHC (as some of the parameter regions where the 
aforementioned signal can be detected overlap with those 
where h\ t 2 — ► a\ai decays might be ineffective in extract- 
ing an NMSSM Higgs signal). 
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